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Abstract 

This work is concerned with multi-dimensional integrals, which are making their appearance in 
few-body atomic and nuclear physics. It is shown that the relevant two- and three-dimensional 
integrals can be reduced to one-dimensional form. This implies that the internal one- and two- 
dimensional integrals can be evaluated in explicit analytic form in term of the familiar generalized 
hypergeometric functions. Some of the integrals are presented here for the first time. 
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I. INTRODUCTION 



One of the effective tools for solving, e.g., iV-body problem both in atomic and nuclear 
physics is introducing the Jacobi vectors (coordinates) which presents a linear combination 
of a standard vectors (i = 1,2,..N) of the particles under consideration. Omitting the 
center of mass motion for a given choice of the Jacobi vectors, the hyperspherical coordinates 
are given by the so-called hyperradius p and by a set fijv-i of angular variables. The latter 
can be expressed through the 2(N — 1) polar angles uji = (#j,0j) of the Jacobi vectors £j 
and (N — 2) hyperspherical angles (pi(i = 2, N — 1). 

In order to estimate the physical properties of the considered system, one needs to cal- 
culate matrix elements of the proper Hamiltonian in the appropriate basis. 

For simplicity let us consider three-particle system (see, e.g., jl, 2|) in the most abundant 
basis set presenting the product 

Rn{p)y{h,l 2 ^},LM ( S l 2), (1) 

where 

m\,m,2 

is the so called hyperspherical harmonic function with definite angular momentum associated 
to the quantum numbers LM. The RHS of Eq.(j2]) contains the function 

7>W% 2 ) = ^^(cos^) Z2 (sin^ 2 )^P; i+1 /^ +1 / 2 (cos2^ 2 ), (3) 

where P^ ,b (z) are the Jacobi polynomials, A/ - ^ 1 '' 2 is normalization constant; the Clebsch- 
Gordan coefficients (Jtimxlzm^LM) , and usual spherical harmonics Yi m (9,(f)). Here \x is 
non- negative integer, and (li,mi) are the quantum numbers corresponding to the angular 
momentum operator associated with the ith Jacobi vector. 

The most commonly encountered form of hyperradial function is the following;: 

R n (p)=C a:n L^(ap)e-^ J (4) 

where Ln (z) are the generalized Laguerre polynomials, and a is a scale factor. The ex- 
plicit forms of the normalization constants C ayU and A/^ 1 '' 2 are not important for the given 
treatment, in contrast to the form of Jacobi vectors 

£l = ^= (2r 3 - ri - r 2 ) , £2 = 7^ ( r 2 - t x ) (5) 



for three particles of equal mass (see, e.g., jl, 3, 4|). Taking the inverse transformation for 
Jacobi vectors (j3J) plus the centre-of-mass vector £ 3 = | (ri + r + r 3 ) and keeping in mind 
that by definition £ 2 = pcos(<^2),£i = psin(<yj 2 ) one obtains: 
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2-r + A v / 3(l-r 2 ) 
2 - r - A v / 3(1 - t 2 ) 



(6) 
(7) 
(8) 



(9) 



where r = cos2<£> 2 , and A presents cosine of the angle between vectors £ x and £ 2 . 
Let us consider the Gaussian-type central potentials of the form 

l> i i>j k 

The examples of such potentials can be presented by the two-nucleon potential models js _ 2|- 
It can be shown that calculations of the matrix elements for the potentials of the form 
fl9]) reduce to evaluating two-dimensional integral of the form 



Z 2 = J (1 - r) /l+1 / 2 (l + r )^+V2p^+i/2 1 / 2+ i/2 (T)p / 1 +i/ 2j / 2+ i/ 2(r)dT x 

poo 

/ LS(x)L^(x)e-^ 2 - x x 5 dx J (10) 
Jo 

and three-dimensional integral of the form 

X 3 = ^ (1 - r 2 )P; i + 1 / 2 ' Zl+1 / 2 (r)P^ +1 / 2 ' i2+1 / 2 (r) C /r x 

/DC ^1 
L^(x)L^(x)e-^- x x 5 dx J e- XKx2 P h (X)P l2 (X)dX, 

where 



11) 



7 = 7l r J 



1 + r 



^ = /3(r) = fl " 



K = KIT) 



(r) = 7rrV3(l-r 2 ), 



2a 



(12) 
(13) 



and Pj (A) are the Legendre polynomials. One should emphasize that the 3-dimensional 
integral (1111) corresponds only to the case of zero angular momentum (L — 0). However, for 
the_potentials considered, the total orbital angular momentum is a good quantum number 
therefore only the wavefunctions with L = were included. 
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It will be shown in the subsequent chapters that both 2- dimensional integral Z2 and 3- 
dimensional integral X3 can be reduced to one-dimensional integrals. The latter implies that 
integrals over A and x variables can be evaluated in the explicit (closed) form. At that, the 
final expressions for the integrals associated to 3-dimensional integral (fTTj) presents a new 
results, at least, from mathematical point of view. 



II. TWO-DIMENSIONAL INTEGRAL 



First, let us consider two-dimensional integral ( ITUj) . Evaluation of the internal integral 
over x can be performed by making use of the known representation for the associated 
Laguerre polynomials 

LS {x)L { £ (x)e-^-*x k dx = J2H C i,l C ll J(k + i + J, 7) (14) 

i=0 j=0 _ 

followed by applying the integral presented on p. 343(2. 3. 15(3)) [8[]: 

J{iy,j) = j\-^ 2 ~ x x v dx = I> + l)(2 1 y^e^D^ 1 f-^j . (15) 

The coefficients introduced in Eq. f[T4l) have a form: 

(n - i + l) k+i (-1)' 

Here, (a) n denotes the Pochhammer symbol, and D a (z) is the parabolic-cylinder function. 
Alternatively, the latter function can be substituted for the Tricomi confluent hypergeometric 
function or Hermite function using relations (13.6.36) or (13.6.38), respectively [9J]: 



In order to derive the latter connection with Hermite function, we additionally made use 
the following relationship (see, e.g., (13.1.29) [9[) 

U(a, 6, z) = z x - h U{a - b + 1, 2 - b, z), (18) 

that will be applied in the next Section too. 

It is worth noting that for integer v = n the parabolic-cylinder function can be presented 
in the form 10| : 



-">--<*> - ^ e ~ 4 £>•>•(")*-' (75) (£) + 

v ^ P2 i(-„-i) e 4 erfc Hn (b\ _ (19) 



where erfc(y) is the complementary error function, m is the binomial coefficient, and i 
represents the imaginary unit. Eq. (fT9l) enables us to derive the additional representation for 
the integral ([Tj 



An ,na (^) 



rzi!rz 2 !(27) 



ni+n 2 +fc 



where 



rS n2 (7) = ^i!n 2 !(2 7 )" 1+ " 2+fe x 



"2 



k+i 



i=0 



z!(A; + z)! 



{n 2 - J + 1) 



3=0 



j!(A; + j)! V 2^7 



i+j+k 



H, 



i+j+k 



2^7 



(20) 



(21) 



n!!n 2 !(27) ni+n2+fc ^ (tix - i + 1, 



v/7 

i+j+k i+j+k 



E 



"2 

E 



(n 2 - j + 1; 



fe+3 



/!(/, . /)! ^ ./!(/■• • ./)! 



x 



£ ( < + { +Ai )(- i )^- + *-. ( 7^=) H-i ( A= ) (22, 



It is worth noting that both T^] n2 and <Sni,n 2 present polynomials (in 7) with the integer 
coefficients, and possess properties: 



lSJ la (0) = 5W na (0) = l. 



(23) 



III. THREE-DIMENSIONAL INTEGRAL 



The aim of this section is to derive the explicit (closed) expression for the internal two- 
dimensional integral 

/oo p 1 

Ll k i \x)L^(x)e-^ 2 - x x k dx j ^ x ™ 2 P h (\)P h (X)d\ (24) 

associated to three-dimensional integral ( ITT]) . To this end, one should first of all to apply the 
well-known Neumann-Adams formula which expresses the product of two Legendre polyno- 
mials as a sum of such polynomials, and then to make use of integral presented by formula 
(2.17.5(2)), p. 428 This yields: 

h 



X ~ 1 \/-T, A hM(- i y i+l2 ~ 2ri h+l 2 -2r+l/2 {kX 2 ) , (h < l 2 ) 



(25) 



where 

(2ji - 2r - l)!!(2r - 1)!!(2Z 2 - 2r - l)!!^ + 1 2 - r)!(2/ 1 + 2l 2 - 4r + 1) 
Wa (/ 1 -r)!r!(/ 2 -r)!(2/ 1 + 2/ 2 -2r + l)!! ' 1 j 

and I m +i/ 2 (-2) are spherical modified Bessel functions of the first kind (m is integer). 

Presenting the product of the Laguerre polynomials in its explicit polynomial form, as it 

was done in the previous section, and making use of Eq. fT2"6]) . one obtains: 




r=0 i=0 j=0 



/■oo 

/ x^' +fc - 1 e-^ 2 -^, 1+Z2 _ 2r+1/2 (kx 2 ) dx, (27) 
Jo 

where coefficients C? k are defined by Eq.( TT6|) . 

Much of what follows are devoted to deriving the explicit analytic expression for evalu- 
ating the integral 

/■oo 

KF{p,n)= I e- px2 - x I^{Kx 2 )x p dx. (28) 
Jo 

We have not found a solution of this problem both in mathematical and physical literature 
even for integer values of p corresponding to Eq. fT27|) . 

Final results that will be obtained here are valid for any real p allows the integral ( 12 8 p 
to be convergent at the given half-integer /i. On the other hand, it is clear that this integral 
converges only for (5 > k. We shall consider real /3 > and k > 0. Note, that the case 
of (3 = k is presented in (see, 2. 15. 6(1), p. 306). However, it is easy to make sure that 
according to definition ( Tl~3|) 

K Ml -r 2 ) . N 

< - = V K 2 _ T < 1 (29) 



for values of — 1 < r < 1. 

First, let us present series expansion for the modified Bessel function of the first kind [9J: 
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2 \ fj, oo i kx 



, x 2k 



Inserting this representation into Eq.f TSHj) . one obtains: 



k=0 



w) ^ r(4fc + v+p + i) / « \ 2 k ( p + i l ii 



(4/3)^^ r(fc + /U +l)fc! V8/3; V 2 '2' 4/3 



Instead of Tricomi confluent hypergeometric function, one can use any of equivalent functions 
as provided by Eqs. (ll5l) - (fT9l . It is worth noting that the cutoff expansion fl3T|) can be 
successfully applied for computing integral (I2"8"j) with any values of (k//3) < 1. However, 
the more value of (ft//3) requires the more length of expansion f[3~Tj) . Therefore, the latter is 
especially effective for very small values of (k//3). 

Let us proceed with derivation of the closed analytic expression for the inte gral (1281) . 
The primary definition for the Tricomi confluent hypergeometric function yields 12] (see, 
7.2.2.(2), p.434): 



V^F 1 F 1 (2k + £±i + //; i; i) ^ ^ (2fc + | + /i + 1; |; i) 



(32) 



T (2^ + 1 + ^ + 1) ^TJ r (2A; + £±i + ^) 

Inserting the first term of the RHS of Eq. (132j) into Eq. (l3~TT) and then changing the order of 
summation, one obtains: 



„ _ ^(^frti + S/i + p+ll/.f 1 F l {2k + ^ + ^ 

^1 - , A n^2+± 1^ 



1. J_ 

' 2' 4/9 



(4/3)^ ^ r(fc + // + l)A;! V8/57 r(2A; + | + /i + l) 

MAL y. 1 ^ r(2fc + ^+/i + n) / K \ 2fc 

2/3^ ^n!r(l/2 + n)(4/3)™^ £;!r(fc + // + 1) \20 

(g) ^ r(n + /i + ^) /n + A* P+l n + p + 3 k 2 \ 

— (2^ — + — — + — (33) 

We used for the latter derivation and we will use in what follows, the well-known dupli- 

n. 

cation formula (6.1.18) |9[ for the gamma functions. 

Next step is making use of the relationship between the Gauss hypergeometric functions 



of the form 2 Fi (a, a + 1/2; c; z 2 ) and the associated Legendre functions of the first kind 
(7.3.1(101), p. 460), which yields for the case of interest: 
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. n + fi p + 1 n + fi p + 3 2 

2^i — ^ — i — — i — r - ; 1 + A*;* 



4 ' 2 4 



2 \ M 



r(M + i)(n) (i-- 2 r^^(^=), (34) 



with z = k//3. 

In this stage, one needs to express the Legendre functions presented in Eq. (l34|) through 
finite sums with the limits (of summation), which are not dependent on n. Note, that 



relationship 3.2(9) 13J enables us to express the associated Legendre functions of the first 
kind through two Gauss hypergeometric functions by different manners, which are presented 
by formulas (3.2.14)-(3.2.31) 13). From the latter representations it is seen that only the 
case of a half-integer \i could give the mentioned above Gauss hypergeometric functions with 
one (of the first two) negative integer parameter, which is not dependent on n. 
Thus, introducing denotation 



jj, — m + 



(35) 



where m is a non-negative integer, and making use representation 3.2(30) 13|, one obtains: 

1 



e 17r U-4, 



V2 



x 



HZ 



1 



2n+p-l 

z \ 4 



l + z 



r(n + p/2 + m + 1) 2 1 V 2' 2z 



2n+p-l 

l + 4 rm + p/2 — m) „ / p l + z 

r(, i+P /2 + i) 2f 'r +i --'" ii+ "n i — 



(36) 



Note, that factor e 1 



" 4 J is correct, but it is different from the corresponding factor pre- 



sented in Ref.|l3|. 

One should emphasize, that due to the factor T(n + p/2 — m) Eq.( l36i) is not valid only 
for the case of even p < 2m (n > 0). This case will be considered later. 

Now, one needs to insert Eq. (I36I) with the hypergeometric functions presented in explicit 
(polynomial) form into Eq. (l34"j) . Then, inserting the result into Eq. fl33|) and changing the 
order of summation, one obtains for the real values of p > —2m — 2, excluding the case of 
even p < 2m: 



vSttk K,\ (-i even p < 2m) 



p/2 



ra=0 



T (m + n + 1 + |) T (n + § - m) 
(2n)!r(l + n+|) (/3 - K ) n 



Fi ( m+ 1, — m; 1 + n + 



p /3 + K 



+ k 



p/2 



00 r (Vi + R ) 



n=0 



(2n)!(^ + «) 



n 



2 2k 

p J3 + k 
2' 2k 



G*-*)-*r(f 



1 r 



p 



m x 



E 



(™+*0!H£r) fc 

- fc!(m-fc)!r(l + fc + §) 



2-^2 



P P. , t 1 P , , , , 1 

2" m '2 +m+1; 2'2 +fc + 1; 4(^) 



(^^(-irE lm+fc) '' (l "^ ( ^ )fc ^ 



k=0 



k\(m — k)\ 



2 '2'4(/3 + k)_ 



(37) 



For the even values of p < 2m, one needs to divide a summation over n by two ranges: 
[0 : m — p/2] and [m — p/2 + 1 : oo], and then - to perform the same procedure, as in the 
previous case. This yields: 

V8ttk /Ci (even p < 2m) = 



/ \ m+l 



P 



r (m + 1) p 



r(| + m + n + l) / n + m+1 p n + m + 2 p 3 k 



n=0 



{2n)\p n 



1F1 



2 
1 



+ 



4' 



(2m + 2-p)\ 



4' m+ 2'/3 2 



+ 



+ K 



m+l 



X 



k=0 



k\ 



3 — p p 
1,1 — k + m;m-\ , m + 2 , 



1 



0-K 



2 

m+l 



2'4(/3 + k) 



+ 



(2m + 1)! x 



k=0 



/c!(m — k)\(m + k + 1) 



3-^3 



1,1, 2 + 2m; 2 + fc + m,m + ^-—^,2 + m- ?; ' 



2'400-k) 



(3* 



One can proceed to consideration of the residual RHS of Eq.( l3TT) . Inserting the second term 
of the RHS of Eq. (1321) into Eq. (1311) and then changing the order of summation, one obtains: 

m+l/2 

/ A r(4fc + 

£^ k\T(k + m + 3/2) \8(3j ^ T (2k + \ + m + l) 



K, 



2m + p + 2) f n\ 2k ^ ( 2k + + "»l I? ^) 



(W 

\ m+l/2 



- r(n + m + ^) 



2r(m + 3/2)/5t+ 1 ^ (2n + l)!/3* 



ra + m p + 3 n + m p + 5 3 k 2 



.(39) 



Once again one should make use of the relationship between the Gauss hypergeometric 
functions of the form 2-^1(0, a + 1/2; c; z 2 ) and the associated Legendre functions of the first 
kind |12| (7.3.1(101), p. 460). However, in this case one sets a = (n + m)/2 + (p + 3) /4, c — 
m + 3/2, z — k//3. Then, one needs to apply representation 3.2(9) [l^ with the parameters 



presented by 3.2(30) 
2F1 



n + m p + 3 n + m p + 5 



131 ] . Thus, one obtains: 
3 
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' 2 






'( 


:::) 








2n+p 


1 - 


z\ 4 


1 + 





; m H — ; z 
4 ' 2' 



T(m+ I 
2^ 



— I— — \ /9\ m+1 2»+p 



l + z^~ T (n+ 2±i -m) 

r(n + ^) 



2 Fi ( m + 1, -m; n + 



T(n + ^ + m) 



Fi (m + l, — m; — n + 



p+3 z+1 
2 ' 2z 

1 — p z + 1 

2 ' 2z 



(40) 
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One should notice that due to the factor T (n+ ^rp- — m) Eg. ( HO"]) is not valid for odd values 
of p < 2m — 1. 

Substituting representation fHUj) with the explicit (polynomial) expressions for the hyper- 
geometric functions into Eq.( |39l) . and then changing the order of summation, one obtains: 



\J 8ttk IC2 f -1 odd p < 2m — 1 



l> {P + K > ^ (m -*)!*! 

fc=0 



K 



2k 



iFi 



p+1 



-k- 



2'4(/9 + «) 



( ^ fi )-^r(^- m )r(^ + «.|x 



(m + A;)!(-l) fc f (3 + k 



" Q [m - k)W (£±3 + k) V 2k 



2^2 



2 J V 2 

p + 1 p + 3 3 p + 3 



m, 



+ m; -, 
2 ' 2' 2 



+ k; 



1 



4(/9 - «) 



For the odd values of p < 2m — 1, one needs to divide a summation over n in Eq. fl39|) by 
two ranges: [0 : m — (p + l)/2] and [m — (p — l)/2 : oo]. And then, it is necessary to perform 
the same procedure, as in the previous case. This yields: 



v8ttk K,2 (odd p < 2m — 1 



m+l 



p+i 



2 vM^J 2 T(m + n+ s ¥) fn + m p+3 n+m p + 5 3/t 2 
— rr 2^ — ^ — "tttttt^^ 2-^1 I — ~ 1 ; — i — ~ 1 ; — ; m H • — 



r(m+|)/3^ ^ (2n + l)!^ 



^ (m + fe)!(^) fc ^ 



4 ' 2 
1 



(2m + 2-p)! 



4 '"" ' 2'/3 2 

m+l 



fc=0 



1,1 — k + m;m -\ , m + 2 



+ 
p 1 



X 



0-K 



m+l 



2'4(/3 + k) 
(2m + 1)! x 



+ 



fc=0 



V 2k / 



fc!(m — k)\{m + /c + 1) 



3-^3 



3 — » » 1 

1, 1,2 + 2m; 2 + k + m,m H -,2 + m - -; 



2'4(/3-«)_ 



A. Final formula for the new integral 



The results obtained above for the integral (128]) with half-integer parameter /x can be 
processed and presented in the compact form. To this end, let us introduce three auxiliary 
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functions: 



Ji(s) = 7i(s,p,m;/3,/c) 



m+l 



r(m+|)/3 



2 



r(m + n+l + £±£) / n + m + l „ +sn + m + 2 p+s 3 k 2 \ 

£ r^TX^ 2i M — 9 — + V> — o — + V- m+ 2 ; ^)' (43) 



ra=0 



(2n + s)\/3 n 



7 2 = 7 2 (p, «) 



1 



(2m + 2-p)! I V + k 



m+l 



fc=0 



fc=0 



jfe! 



3 — p p 
1, 1 — k + m; m H — , m + 2 — - 



2'4(/3 + k) 



(2m + 1)! 

(/3 - k)™* 1 



x 



fc!(m — k)\{m + + 1) 



3^3 



3 — p pi 
1, 1, 2 + 2m: 2 + fc + m, + m, 2 + m - -; 



2'4(/3-k) 



-1^ 



^ (m + fc)!r(^-fc) ^ + if . 



7 3 (s) =jF3{s,p,m;p,K) = 
p + s 



k=0 



(m — k)\k\ 



2k 



2s + 1 



4(/3 + k) 



+ 



r (2±2 - m) T (E±2 + m+l) ™ 



fc=0 



m 



(m + A;)!(-l) fc + « 

jfc)!fc!r(2p + fc + l) V 2k 



2-^2 



P + s P + s 2s + lp + s 1 
m, — h m + 1; — - — , — h k + 1; 



2 " ' 2 ' "" ' "*"' 2 ' 2 ' ' *'409-k)_ 
The integral of interest can be expressed in term of these functions as follows: 

7i (0) + 7*2 - 7*3(1) even p < 2m 
-7i(l) -7 2 + 7 3 (0) oddp<2m 
73(0) — 73(1) otherwise 



I) ' ' - V87TK 



(46) 



It is seen that representation ( )46l) cannot be used for very small values of parameter k. How- 
ever, in this case, the cutoff representation ( 13TT) can be applied with advantage. Moreover, 
one should notice that the closed form ( 14"6"|) presents the difference of two large quantities, 
because /Ci and /C 2 have different signs, but their absolute values are very close. At that, the 
relation |/Ci/(/Ci — /C 2 )| increases very quickly with parameter m of the spherical modified 
Bessel function of the first kind. Therefore, the quantities fC\ and /C 2 (i.e., functions 7j) 
have to be calculated with high accuracy, especially for large m. 
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It is worth noting that a few-body problem is related to the integral (jlSl) with integer 
power p and even m, whereas formulas P3]) - (H6]) are valid for any integer m > and any 
real p > —2m — 2. 
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